Scattering theory of magnetic/superconducting junctions with spin active interfaces 
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We formulate a generalized scattering field theory a la Biittiker describing particles transport 
in magnetic/superconducting heterostructures. The proposed formalism, characterized by a four- 
component spinorial wavefunction of the Bogoliubov de Gennes theory, allows to describe the spin 
flipping phenomena induced by noncoUinear magnetizations in the scattering region. As a specific 
application of the theory, we analyze the conductance, the magnetoresistance and the generation 
of spin-torque produced by an applied voltage in a spin-valve system. Quantum size effects and 
quantum beating patterns both in the conductance and in the spin-torque are carefully described. 



I. INTRODUCTION 



Nanoscale structures involving Normal (N), Ferromag- 
net (F) and Superconductor (S) junctions, the so-called 
heterostructures, involve interplay of superconducting 
and ferromagnetic order parameters providing a novel 
opportunity to study the influence of the spin degree of 
freedom on transport and thermodynamic properties of 
such systems. A paradigmatic example is represented 
by the normal metal (N) /superconducting (S) bilayer. 
In this system the sub-gap transmission of an electron 
propagating from the N-side towards the S-side is for- 
bidden due to the absence of available electronic states 
within the superconducting gap. Thus, in order to con- 
serve the charge current, a propagating Cooper pair is 
transmitted in the superconductor, while an hole is re- 
flected in the normal metal. This anomalous reflection, 
described by the pioneering work of Andreev^ , has been 
successively described by Blonder et al^ in the language 
of scattering theory (the so-called BTK theory) . Differ- 
ently from methods based on the transfer Hamiltonian 
formalism (THF), the BTK scattering approach does not 
assume weak coupling approximation and allows to study 
transport in NS heterostructures from the metallic (high 
transparency of the interface) up to the tunneling limit 
(low transparency). In the tunneling limit, the agreement 
between BTK theory and the Green's function approach 
shows that, in the absence of many-body correlations, a 
scattering theory is extremely suitable for studying trans- 
port in heterostructures avoiding time-consuming meth- 
ods. Almost ten years after the BTK formulation, a de- 
velopment for F/S interfaces based on a BTK-like theory 
able to describe the Andreev reflection physics within 
the scattering approach appeared^. A further improve- 
ment of the original BTK formalism has been succes- 
sively brought out by Anantram et alA who reformulated 
the BTK theory in terms of the scattering fleld theory^ 
(BSFT) originally conceived by Biittiker^ for mesoscopic 
(normal) systems. The formalism provides a clear cor- 
respondence between physical observables and the scat- 
tering matrix of the system within a second quantization 
formalism, avoiding the by-hand construction character- 
izing some parts of the BTK theory^. The Anantram 



and Datta work^ is based on a two-component spinorial 
Bogoliubov-de Gennes^ theory (as in the original BTK 
theory) which is a convenient representation in the anal- 
ysis of spin conserving processes. However, recently, the 
need to explore the interplay between superconductivity 
and magnetism stimulated the realization of magnetic su- 
perconducting heterostructures. In the simplest case of 
a ferromagnetic F/S bilayer^, the bulk magnetization of 
the F-side can differ from the one at the interface and spin 
flipping phenomena can take place at the F/S junction. 
In order to fully describe the physics of the magnetic su- 
perconducting heterostructures one needs to generalize 
the BdG formalism to a four-component spinorial repre- 
sentation. To the best of our knowledge, a scattering field 
theory a la Biittiker for magnetic superconducting het- 
erostructures which properly treat spin active phenom- 
ena at the interface is not yet available and our work aims 
to fill such vacancy. In particular, the second quantized 
form of the BTK theory offers important advantages in 
obtaining non-local quantum properties (i.e. correlation 
functions) which are not provided by the BTK theory. 
Following the work of Anantram et alA, the construction 
of a spinful theory is important to correctly describe in- 
terplay phenomena between superconductivity and mag- 
netism. 

For the presentation we follow the formalism of 
Refs.0|j] and focus our attention on quasi-one- 
dimensional systems, while the generalization to the 2- 
dimensional case is left for a forthcoming work. Within 
our generalized scattering formalism we derive the ex- 
pression for the charge and spin currents and the respec- 
tive linear response to an external bias, i.e. the conduc- 
tance and the spin-torque. The second part of the paper 
is devoted to the study of the spin polarized transport 
in a spin-valve system by the developed formalism. In 
particular, we focus on the evidence of the Andreev re- 
flection in the subgap transport at varying the magnetic 
interaction and on the quantum size effects due to the 
interlayer width to make a link to experimental obser- 
vations. Among the transport properties, we do analyze 
the conductance, the magnetoresistance and finally the 
spin torque as a probe of spin-polarized transport. 

The organization of the paper is the following: In SecHIl 




Scattering region 



FIG. 1: Magnetic superconducting heterostructure as de- 
scribed in the main text. The sequence of materials 
[Qi|a2|...|an] can be a„ € {F,S,NM}. 



we introduce the Bogoliubov-de Gennes Hamiltonian and 
present the scattering field theory generalized to include 
spin flipping phenomena. We then derive the expression 
of the charge and spin current and of the respective differ- 
ential conductance in the presence of an external bias. In 
Sec lIIII we present the results of the linear response ob- 
servables: the conductance, the magnetoresistance and 
the torkance for the structure shown in Figl2] Conclu- 
sions and perspectives are given in Sec. IV. 



II. 



THE MODEL AND FORMALISM 



We consider a one-dimensional magnetic super- 
conducting heterostructure connected to normal 
nonmagnetic leads (see FigllJ. The scattering region 
is made by magnetic and superconducting regions a^ 
arbitrarily disposed along the transport direction. The 
system is conveniently described by a 4-component 
Bogoliubov-de Gennes (BdG) formalism in which 
the quantum state of the system is described by the 
wave- function |*(x,t)) = X^^.a <^(9CT(a;,i)|/3) ® \<t), 
where /3 S {e, h} is the particle index and a € {+, —} 
represents the spin orientation along the quanti- 
zation axis. Introducing the standard notation of 
the BdG theory we set 4>ea(x,t) — >■ Ua-{x,t) and 
4>ha{x,t) — )■ Va{x,t) thus the state vector can be put in 
the form |^(a;,i)) ~ {u^{x,t),ui{x,t),v^{x,t),vi{x,t)y . 

Using a tensor product notation (see APPENDIX \^ , 
the quasiparticle Hamiltonian can be represented as^ 

H = P,,® (H, + U{x)) + Phh ® [-{H, + Uix))*] 

+ Peh ® A{x) + Phe ® Hx)\ (1) 

where P^p = \ct){P\ (a,/3 = e,h) are particle and/or 
hole projectors with |e) = (1,0)* and \h) = (0,1)*; He, 
U{x), A(x) are operators written in the spin-space basis 
\a = ±), \+) = (1,0)* and |-) = (0,1)* (i stands for 
the transposed vector) . Specifically, He is kinetic energy 
operator defined by 

^'^' z. 1 



He 



2m 



Is 



(2) 



being Ep the Fermi energy and I^p the identity operator 
in the spin space. The potential energy U{x) may include 



spin dependent Zeeman or spin-orbit coupling terms (e.g. 
U{x) = h{x)-a) and finally for an s-wave superconductor, 
A(a;) — i(7yA{x) where ay is a Pauli matrix. In the 
following we characterize the superconducting regions by 
a constant order parameter A (i.e. we assume a step-like 
behavior of the gap at the N/S interfaceaiS) . 

Within our tensor product representation, the charge 
and spin current operator in first quantization can be 
written in a compact form as: 



J,- 



2m 



d. 



^. 






(3) 



where the index /i = 0,1,.. .,3, so that Jq is the charge 
current operator, while Ji, J2 and J3 represent the three 
space-components of the spin current operator. In the 
above notation wc introduced go = 9e, 91 = 92 = 93 = 



1 and (a'}, 



ct:; 



cr,; 



CT, 



/3 



ripls 



r]h — ± aiiu yu^) — u^ — Up, Uq — iip'^sp, 

|e| the electron charge and a^ the ^-th Pauli 



h/2, 7]e - 

being qe 

matrix. 

The above expression for the currents derives from the 

conservation laws of the charge Q and spin 5"^ densities. 



Q = Qe^VpP^ 



PP 'i9 lisp 



(4) 



Sf, = (?i/2)^P;3/^®a/^ 



In particular, the charge and spin densities conserva- 
tion law can be put in the form of a continuity equation 
with source/sink terms. Concerning the charge continu- 
ity equation, the source/sink term is related to the di- 
vergence of the Cooper pairs current^!, while in the spin 
density case such term is related to the spin-torque. In- 
deed, a spin-torque T^ can be generated by a magnetic 
potential of the form U{x) — h{x)-a and its /i-component 
is represented by the operator 



T„ 



Y^ppp® 



h{x) X 'a^ 



(5) 



The continuity equation for the spin density can thus be 
written as dtSu + dxJu = Tu- 



A. Scattering formalism 

In the following we describe the scattering field theory 
for magnetic/superconducting heterostructures. In con- 
structing the theory, the Andreev approximation (which 
neglects the difference between the particle and hole mo- 
mentum) is performed only in the external leads, while 
inside the scattering region the exact wave-functions are 
considered. This approach is performed to properly treat 
the phase-coherent phenomena in the scattering region 
and to correctly capture all the Andreev-reflections prob- 
abilities. 



The scattering field (see APPENDIX |B]) in the j-th lead 
can be written as^^: 



^jix,t) 



dE exp[—iEt] 
y/2Tihv{E) 



E 

a'j0iE)e'''^^ + b' 



(6) 
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iE)e- 



-ikffXi 



where kp — ripk{E), while the scattering operator a'^-JE) 
{Mp{E)) destroys an incoming (outgoing) particle of 
species /3 S {e, h} and spin projection a G {+, — } 
in the lead j and within the Andreev approximation 
Vji3a{E) « v{E) ^ hk{E)/m. The scattering field de- 
fined above generalizes the one introduced in Ref.[4| to 
the spinful case. The outgoing field operators are re- 
lated to the incoming field operators by the scattering 
matrijiHiii: 



b-sit) 



E 

i'a'13' 



s: 



PP' 



,{M'p'it)- 



(7) 



Since the scattering states given by ([5]) form a complete 
set of mutually orthogonal states (completeness relation), 
the fields 6f satisfy the canonical commutation relations 

{^i/3' (^^'fl')^} = 5ii'Sj3pi5aa' and the current conservation 
ensures the unitary condition of the 5-niatrix: 



i'cr'P 



qPP ( qD> 



bli' 



SikSbpSsa- 



(8) 



The quantum statistical properties of the leads are 
defined by the expectation value (aJj^(-E)af^(£")) = 
S^J6sJapS{E - E')fja{E), fja{E) being the Fermi dis- 
tribution of the particle of species a in the electrode j. 
The field operator ^j acts on a many-particle state and 
thus the expectation value Oj of the generic operator O 
in the j-th electrode is given by Oj = {'^.O'iij), where 
the notation (• • •) stands for the quantum statistical av- 
erage. 



B. Charge current and differential conductance 

In this Section we derive the two-terminal conduc- 
tance of the magnetic/superconducting heterostructure 
depicted in FigHJ For a multi-terminal device, the av- 
erage charge current Jq flowing through the i-th lead is 
given by Jg = (^J(a:,i) Jo^i(a;,i)) and using the field 
representation ([6]), it can be expressed in terms of the 
scattering matrix as 



^ii = fx+ {-fqeV/2, being ^ = {fii + /X2)/2 and V the 
bias voltage. Taking as zero of the energies the electro- 
chemical potential of the scattering region /is, we can 
write fj^{E) = f{[E + Tj^^ifi, - /i,)]/(KBT)), being T 
the temperature. Let us note that fig — (p-i + A'2)/2 
only in the symmetric case. In the nonsymmetric case 
/is 7^ (/ii + M2)/2 — /i and thus /is must be deter- 
mined self-consistently to conserve the charge current 
(i.e. J2i "^o(^'A*s(V')) = 0) as described in the Appendix 
[Cl Within the linear response theory, the charge cur- 
rent flowing through the i-th lead is obtained as li = 
J2j 9ij{^''j ~ ^^s) — GiV where gij is the conductance ten- 
sor whose expression (see Appendix [C|) is: 



gik 



e 






(10) 



.Pa, 



where the sum rule J^ja-^ii^i-^) ~ -^ ^^^ been used. 
For a generic structure the two-terminal conductance in 
terms of the conductance tensor is given by: 



G 



.922.911 — 921912 



Z^ij 9ij 



(11) 



Let us note that in the symmetric case the above relation 
can be simplified as Gsym — [gii ~ .9i2)/2 and thus: 



G 



sym 



^^Jd^i-ddiOUx 



(12) 



Mim + M'l^iO + M'ltiO + Mf,m- 



C. Spin currents and spin-torque 

When the heterostructure contains magnetic layers 
apart the superconducting ones, the application of an 
external bias produces a spin current coexisting with the 
charge flux. Differently from the charge, the spin density 
is not conserved and thus a spin torque is exerted along 
the nanostructure. In particular the spin torque results 
from the divergence of the spin current as discussed in 
Ref.[l3. 

In order to derive an expression of the spin-torque, 
we need first to derive an expression of the spin cur- 
rent in terms of the scattering matrix. This is given by 
the quantum average of the spin current density operator 
J* = (^P] J^^&i) (/i = 1,2,3) and using (6) one obtains: 



paj 



dE 



iPa 



2S,,Sa.p-M':"iE) f.^iE), (9) 



where 7Wf/(^) 



Tr[S^fiE)Sl''m, Tr[ 



indicates 

the trace on the spin indices, while S^°'{E) are matri- 
ces with respect to the spin indices. When a symmet- 
ric potential drop is applied to the system, the elec- 
trochemical potential in the i-th lead can be written as 



Jl 



E 

aPj 



dE 

47r 



■^Pa, 



Tr[S^;\E)a^^S':;{E)]f,^{E), (13) 



where J' represents the /i-component of the spin current 
density generated in the «-th lead. Let us note that 
Ea. (IT5)) represents an expression of the spin current 
beyond the linear response regime. Moreover in calcu- 
lating the spin current, the charge current conservation 



through the system must be monitored. In fact the 
spin current is not conserved due to the presence of 
a spin-transfer torque acting on the local magnetic 
momentum of the magnetic region and thus a violation 
of the charge conservation law may artificially change 
the spin current gradient. 

In order to relate the spin current gradient to the spin 
torque we can consider the continuity equation of the 
spin density S'^: 



dtSn 



Jfi — -L^- 



(14) 



Under stationary condition, i.e. dtS^ — 0, one can apply 
the Gauss-Green theorem to Eq. ([T4|) . Let us consider a 
cylindrical surface of volume fl encircling the scattering 
center and with axis collinear to the transport direction 
x: 



V • Jf,dV 



T^dV 



Jfj^x ■ ds. 



(15) 



where E is the cylindrical surface S — si U S2 U s; with 
si and S2 the areas collinear to the transport direction 
(si = —§2 — x) and s; the lateral surface of the cylinder. 
Since for a quasi-one dimensional system, the physical 
quantities J^ and T^ can be considered uniform along 
the radial direction, one obtains: 



E 






SR 



Tf^dx = 0, 



(16) 



where the integral is performed over the scattering region 
(SR). Eg. dTBl) is the Kirchhoff's law for the spin current 
and it can be used together with Eg. p^ to derive the 
total spin torque t^ = J^^ T^dx produced along the sys- 
tem: 
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dE 

47r 



Tr[Sl-\Ey 



lia 
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fJ- 1-3 



{E)]f,^{E). (17) 



However, in spin valve devices it is important to resolve 
the spatial dependence of the spin torque density. Indeed, 
a spin valve is a device made by two magnetic layers sepa- 
rated by an interstitial nonmagnetic region (spacer). The 
relative volume of the magnetic layers can be taken very 
different: one is the pinned magnetic layer (fixed layer), 
having the largest volume, and the other is the thin mag- 
netic layer called free-layer (FL). When a spin-polarized 
current interacts with the thin ferromagnetic layer it un- 
dergoes a spin-filtering and the result is, in general, that a 
spin-transfer torque is applied to the magnetic layer. The 
energy required to change the magnetization direction of 
the free-layer can be provided by the flux of spin polar- 
ized current activated by the application of an external 
driving field (e.g. dc voltage bias or ac modulations). 

In order make an explicit calculation we consider a spin 
valve, as the one of Fig. 2, with a Zeeman potential of the 
form: 



U{x) = [^5{x)hi + h{x)n2] 



(18) 
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FIG. 2: Superconducting spin valve as described in the main 
text. 



where fii — {sm{6) , , cos{9)) represents the magnetiza- 
tion direction of the FL, n2 = (0,0,1) is the direction 
of the magnetization of the fixed layer, while h{x) is a 
step-like function. In this model the free-layer is directly 
connected to the first lead (i.e. the region x < 0) and thus 
the spin torque acting on this region is simply related to 
the non-equilibrium spin density produced by the bias in 
the first lead. Using Eq.®, the zero temperature spin 
torque components acting on the free-layer in the linear 
response regime can be written as follows: 



Til 



Ti = 



eVT 



Att 



J2Tr[S^,r''^'yS?r]^o.X, 



(19) 
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Y,Tr[S',f {sm{0)<j^^ - cos(0)af )5^"]^aA,, 



Q/3j 



where F = {kp'y)/Ep, while the coefficients Xj are given 
in APPENDIX O Furthermore, T\\ and T± are the com- 
ponents of the spin torque parallel and perpendicular to 
the plane of the magnetization of the fixed layer (i.e. 
T = T\\v\\ + T±u±) whose directions are defined by the 
vectorslS,: 



"W 

V] 



'XC0s{6) + zsin(6') 



(20) 



III. 



SPIN-POLARIZED TRANSPORT IN 
SPIN- VALVE SYSTEMS 



The scattering formalism developed so far can be em- 
ployed to describe the linear response properties (i.e. 
conductance and torkance) of the superconducting spin 
valve depicted in FigH (sec also APPENDIX |D]). The 
system is described by the BdG Hamiltonian given 
in Eq.(IT]), where the superconducting gap operator is 
A(a;) = iayA9{x)9{di — x), 9{x) being the Heaviside 
step function while the two magnetic regions, namely Fl 
and F2, are modeled by the Zeeman potential given in 
Eq.dlHl), where h{x) = Eph^Oix - di)9{di + d2 - x). 
Furthermore an additional barrier potential of the form 
Usix) = J2j=0,l,2^3^i^ ~ ^j)^'>p (xo = 0, cci = di. 



X2 = di -\- d2) is introduced at the interfaces. The 
Vj are related to the dimensionless BTK parameters 
Zj — 2raVjl{h^kF) which measure the interface trans- 
parencies, and are given by the transmission and reflec- 
tion probabilities 7j, TZj via the relation Zj = \jT^j /Tj ■ 
Finally, the s-wave order parameter is taken in dimen- 
sionless form as ?7 = A/Ep and its relation to the BCS 
coherence length ^ is given by kp£^ = 1/ry, ^ sa hvp/i^A). 
In the following we set 77 = 1/200 and kp ~ lA~^ which 
are suitable phenomenological values for conventional su- 
perconducting materials such as Nb- ^. 

In the subsequent analysis, adopting the same line of 
Ref. [J| , the self-consistent computation of the supercon- 
ducting order parameter is neglected. In fact we will con- 
sider the low-bias regime (i.e. eV/A <^ 1) under which 
the spin accumulation in the superconducting region is 
unable to produce a relevant suppression of the super- 
conducting gap. 

A different mechanism of modification of the super- 
conducting gap could be induced by the size of the su- 
perconducting region as reported in Ref. [10]. However, 
as shown in Fig. 2 of that work, the superconducting gap 
saturates to the bulk value as a function of the thick- 
ness of the superconducting layer already at values of 2^ 
(^ being the BCS coherence length). These features are 
quite generic and seem to be robust for any value of the 
scattering potential at the F/S interface and also for par- 
allel or anti-parallel magnetizations in the ferromagnetic 
leads. Thus we conclude that in our analysis neglecting 
the self-consistency of the gap does not induce quantita- 
tive important changes. 



A. Differential conductance and magnetoresistance 

In Fig|3] we report the differential conductance G as 
a function of the energy e/A, with e = eV, computed 
setting the model parameters as follows: kpdi = 600, 
kpd2 = 250, zi = Z2 = Z3 ^ 0.1, 6* = 0, 77 = 1/200, 
hz = —0.45. At increasing the Zeeman interaction F of 
the thin layer, a lowering of the conductance is observed 
below the gap. This is due to the fact that Andreev 
reflection processes dominating the transport properties 
below the superconducting gap become suppressed. 
This behavior qualitatively reproduces the experimental 
observations reported in Ref. [16] by STM technique. 
The effect of the spin active barrier on the transport 
properties of the system is analyzed in FigsH) The 
figures represent the differential conductance G com- 
puted using the parameters: kpdi = 575, kpd2 = 350, 
zi = Z2 = Z3 = 0.1, r = 0.85, ri = 1/200. For both 
the upper and lower panel the full line is computed by 
setting — 7r/2, while the dashed line is obtained fixing 
= 0. The Zeeman energy of the fixed layer is taken 
hz = —0.5 in the upper panel and h^ = —0.75 in the 
lower one. The analysis of the figures shows that the sub 
gap transport is not much sensitive to the magnetization 
direction, while the quasi-particles transport depends 




FIG. 3: Two-probe differential conductance G (in unit of 
e^//i) as a function of e/A for the model parameters; kpdi = 
600, kFd2 = 250, 21 = 22 = 23 = 0.1, S = 0, 77 = 1/200, 
hz — —0.45. The parameter F takes values ranging from 0.9 
(top curve) up to 2.2 (bottom curve) and it is increased with 
constant step of 0.1 going from top to bottom curve. 



on the orientation of the magnetization of the free-layer 
and more harmonics appear in the oscillating behavior 
of G above the gap. We do also observe a lowering of 
the differential conductance as a function of hz from the 
upper to the lower panel. The origin of the oscillations 
above e ~ A is due to releasing the Andreev approxima- 
tion and are related to the formation of quasiparticles 
resonances above the gapv^i. 

In order to describe the magneto-transport properties 
of the system, we introduce the magnetoresistance 
(MR) defined as follows: MR = [Gp - Gap]/Gap, 
where we defined Gp = G{\hz\,T,9 — 0) and 
Gap = G{-\hz\,T,e = 0). In Fig [S] we report the 
MR as a function of e/A setting the remaining parame- 
ters as follows: kpd2 — 250, zi = Z2 = Z3 = 0, F = 0.5, 
77 = 1/200, e = 0, \hz\ = 0.5. The different curves are 
related to different width of the superconducting region 
and in particular the dashed line indicates kpdi = 600, 
the full line kpdi — 800, while the dashed-dotted line 
kpdi — 1000. The analysis of the results shows a change 
of sign of MR as a function of kpdi for e < 0.5A which 
indicates a change in the relative magnitudes of Gp 
and Gap. Furthermore, the small subgap values of 
the MR indicate the inefficiency of the spin polarized 
transport operated by the Cooper pairs. On the other 
side, above the superconducting gap the quasi-particles 
transport efficiently provides spin polarized currents 
and thus MR values ranging from —10% up to 20% 
are observed. The low values of MR below the gap 
are due to a thickness of the superconducting region 
larger than the coherence length ^. Indeed, the curves 
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FIG. 4: Two-probe differential conductance G (in unit of 
e^ /h) as a function of e/A for the model parameters: kpdi — 
575, kpdz = 350, zi ^ Z2 = zs = 0.1, T = 0.85, r; = 1/200. 
The full line is computed by setting 6 = 7r/2, while the dashed 
line is obtained fixing = 0. The Zeeman energy term of 
the fixed layer is taken hz = —0.5 in the upper panel and 
hz = —0.75 in the lower one. 



in FiglS] are obtained for di > 3^, i.e. for a thickness 
such that the quasi-particles current coming from the 
normal leads is almost fully converted in non-polarized 
supercurrent. The latter point is evident in FiglS] which 
presents the MR as a function of the size kpdi of the 
superconducting region computed for the following set 
of parameters: e/A — 0.01, zi — Z2 ^ z^ — 0, T = 0.5, 
77 — 1/200, \hz\ — 0.5, 9 = 0. As the superconducting 
spacer becomes larger than 3^ (i.e. kpdi = 600) a strong 
suppression of the MR is observed for all the curves, 
while below this threshold the quasi-particles current 
is not efhciently converted in unpolarized supercurrent 
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FIG. 5: Magnetoresistance MR as a function of e/A for the 
model parameters as follows: kFd,2 = 250, zi = Z2 = zs = 0, 
r = 0.5, ?7 = 1/200, e = 0, \hz\ = 0.5, while kpdi = 600 
(dashed line), kpdi = 800 (full line) or kpdi = 1000 (dashed- 
dotted line). The inset contains the MR behavior within the 
energy range [0,0.5]. 



leading to a residual polarization responsible for sizeable 
values of MR (w ±10%). The latter resuhs imply that 
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FIG. 6: Magnetoresistance MR as a function of the size kpdi 
of the superconducting region computed for the model pa- 
rameters: e/A = 0.01, zi = Z2 = Z3 = 0, T = 0.5, 77 = 1/200, 
\hz\ = 0.5, 9 = 0, while kFd2 = 250 (circle, o), kFd2 = 300 
(square, D) or kFd2 = 400 (diamond, 0). The inset contains 
the MR vs kpdi in the range [850, 1200]. The sampling step 
is 25 (i.e. 2.5nm). 

a competition between superconducting and magnetic 
properties becomes relevant for di < 2^, i.e. di < 40nm 
for Nb superconductors. This result is consistent with 
that found in Ref. ,5,] and the above conditions can be 



easily tackled in nanostructured devicesii. 
Regarding the oscillations observed they come from the 
formation of resonant states below the gap and their 
period is of the order of the coherence length. From the 
analysis above it is evident that the behavior of the MR 
is related to the amount of polarized current transmitted 
to the free-layer. This quantity on turn depends on (i) 
the efficiency of the fixed magnetic layer in polarizing 
the particle current and (ii) on the transmission of 
the polarized current produced by the polarizer (i.e. 
the fixed layer) through the spacer region. Point (i) 
is investigated in FiglT] where the MR is reported as a 
function of the size kpdi of the superconducting region 
for different values of \hz\ setting the model parameters 
as follows: e/A = 0.01, zi = Z2 = Z3 = 0, T = 0.5, 
rj = 1/200, 9 — Q, fci?(i2 = 650. Apart from the general 
aspect similar to the one of Fig|6l one observes that an 
increasing of the Zeeman energy hz of the fixed layer 
produces higher values of MR for a superconducting 
spacer width smaller than kpdi — 400. Furthermore, 
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FIG. 7: Magnetoresistance MR as a function of the size kpdi 
of the superconducting region computed for the model pa- 
rameters: e/A = 0.01, 21 = Z2 = 23 = 0, r = 0.5, r? = 1/200, 
61 = 0, kpdi = 650, while \h^\ = 0.1 (circle, o), \h^\ = 0.25 
(square, D) or \hz\ — 0.5 (diamond, 0)- The inset contains 
the MR vs kpdi in the range [850, 1200]. The sampling step 
is 25 (i.e. 2.5nm). 



increases at increasing the barrier heights Zj from up 
to 0.1. In higher dimension (2D or 3D) we do expect 
that interface disorder can reduce the amplitude of 
such oscillations. Finally the quantum size effects are 
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FIG. 8: Magnetoresistance MR as a function of \hz\ computed 
for the model parameters: e/A = 0.01, z\ — Z2 = z^ = 0, 
F = 0.5, r? = 1/200, 61 = 0, kFd2 = 650, kpdi = 200. 

displayed in Fig l9] where a density-plot of the MR in the 
plane {kpdi,kpd2) is shown for the model parameters: 
e/A = 0.01, zi = Z2 = Z3 = 0, r = 0.5, ry = 1/200, 
9 — Q and \hz\ = 0.3 (upper panel) or j/i^j — 0.5 (lower 
panel). The dimensionless sampling step adopted in the 
numerical simulations is 10 (i.e. 1 nm) for the upper 
panel and 25 (i.e. 2.5 nm) for the lower paneli^. The 
overall behavior of the curves presented in FiglH] show 
oscillating patterns and a change of sign of the MR as 
a function of the geometric parameters of the system. 
The comparison between the upper and lower panel 
shows the effect of the magnetic energy \hz\ of the fixed 
layer in rotating the wave front of the curves. This is 
particularly evident for the MR as a function of ^^^2 
(i.e. the length of the fixed layer) for a fixed size of the 
superconducting spacer. This dependence can represent 
a relevant information for the experiments. 



the behavior of the MR as a function of hz is expected 
to be proportional to Thz cos(6'), i.e. the scalar product 
of the magnetic momenta of the ferromagnets. This 
is found in FiglS] where we plot the MR as a function 
of \hz\ for the other model parameters: e/A = 0.01, 
zi = Z2 = Z3 = 0, F = 0.5, 77 = 1/200, 6* = 0, 
kpd^ = 650, kpdi — 200. The analysis of the figure 
shows a linear behavior with respect to \hz\ with a slope 
proportional to rcos(0), while an additional oscillating 
pattern is observed. Such superimposed oscillations 
depend on the interface potentials and their amplitude 



B. Spin-torque 

Up to now we focussed our attention on the MR; how- 
ever an additional probe of the spin polarized transport 
through the system is provided by the spin torque (T\ 11) 
acting on the free-layer, see Eg. p^ . Despite the few ex- 
perimental reports concerning the direct measurement of 
this observable, recently difficulties in making quantita- 
tive measurements of the spin-torque seem to be over- 
come. In particular magnitude and direction of the spin 
torque have recently been measured in magnetic tunnel 
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FIG. 9: Magnetoresistance MR as a function of kpdi and 
kFd2 for the model parameters: e/A = 0.01, zi = 22 = 23 = 
0, r = 0.5, T) = 1/200, e = and \K\ = 0.3 (upper panel) 
or I fez I = 0.5 (lower panel). The density- plot shows the large 
scale structure of the oscillations. 



junctioniS' leading to a substantial understanding of the 
angular momentum transfer in these systems. These de- 
vices are of primary interest for the applications and rep- 
resent excellent probes of the possibility to electrically 
control (using dc or ac signals) the magnetic degrees of 
freedom (i.e. the free layer magnetization). Within this 
framework, the study of superconducting spin valves (as 
the one depicted in Fig. 2) can clarify the mechanism in- 
volved in the angular momentum transfer through a thin 
superconducting layer thus constituting a complemen- 
tary tool in investigating the interplay between super- 
conductivity and magnetism. A systematic analysis of 
these structures, also including different symmetries of 
the superconducting order parameters, could be useful 
to probe exotic pairings and their ability in supporting 
spin polarized current. 

To start this study, we perform an analysis of the s- 
wave case here. In the following we take eV — >■ 1 and 
thus the quantities 2\.|| in units of eT^ coincide with the 
derivative of the spin torque with respect to the bias in 
the linear response regime, i.e. the so-called torkance. In 
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FIG. 10: Parallel (dashed line) and perpendicular (full line) 
component of the spin torque Tx^n as a function of 6* computed 
for the model parameters: e/A = 0.01, h^ = 0.5, F = 0.5, 
77 = 1/200, e = 7r/2, kpdx = 300, ^^^2 = 300, while z = Q.l 
for the upper panel, 2 = 0.3 for the middle panel and z = 0.6 
for the lower panel. 



FigHniwe report the parallel (dashed line) and perpen- 
dicular (full line) component of the spin torque 2^ 11 as 
a function of 9 computed setting the model parameters 
as follows: e/A = 0.01, /i^ = 0.5, F = 0.5, 77 = 1/200, 
9 = 7r/2, Rpdi = 300, kpd2 = 300, where we use z — 0.1 
for the upper panel, z = 0.3 for the middle panel and 
z = 0.6 for the lower panel. The spin torque components 
present an almost sinusoidal behavior as a function of the 



magnetizations angle 9 and thus exhibit vanishing values 
for 9 = 0,7r. The maximum values of Tj_\i are observed 
close to 6* = ±7r/2. By analyzing FigfTOlwe do observe an 
increasing of the maximum value of T± and a change of 
sign of T\ I going from the upper to the lower panel (i.e. by 
increasing z from 0.1 up to 0.6) . The latter behavior is 
attributed to the difference of spin polarized currents at 
the interface. The maximum (minimum) value of Tj_ in 
the lower panel (see Fig fTO)) close to ^ = — 7r/2 {9 — 7r/2) 
takes an absolute value of O.lfieV in the presence of an 
applied bias of 1.5meV. This value is of the same order of 
magnitude of that obtained in the case of nonsupcrcon- 
ducting spin- valves (see for instance Ref.^]). This fact 
points out that nanostructured superconducting material 
can support a spin polarized particles transport in agree- 



ment with recent experimental findings, see e.g. Ref. 17 1. 
However the values of the spin-torque strongly depend on 
the interface properties, i.e. on the parameters Zj in our 
model, and thus a comparison with the experimental data 
can be done only by considering Zj as phenomenological 
fitting parameters. The behavior of the spin torque 2\ ii 
as a function of z (where we set zi = Z2 = z-^ = z) for 
different thickness of the SC layer is shown in FigfTTlfor 
the choice of parameters: e/A = 0.01, hz = 0.5, F = 0.5, 
77 = 1/200, 9 = tt/2, kFd2 = 300. All the curves present 
maximum values of the torkance close to z ~ 1, while the 
maximum value of the spin torque component T± is in the 
range 1 — 2/ieV. For highest values of z, i.e. z > I, Tlji 
start to decrease as an effect of vanishing particles flux 
through the interfaces. On the experimental side, all our 
analysis can represent an efficient way of detecting the 
spin polarized effects in the magnetic/superconducting 
heterostructures, despite the experimental difficulties of 
engineering reproducible interfaces. 



IV. CONCLUSIONS 

In this work a scattering field theory for quasi-one- 
dimensional magnetic heterostructures containing s-wave 
superconducting regions has been developed. The second 
quantized form of the scattering fields for the spinful case 
allows a direct link between physical observable and their 
relation to the scattering matrix describing the system. 
Our formalism take fully into account Andreev reflec- 
tions in the presence of spin-flip phenomena. We for- 
mally derived spin and charge current and all the quan- 
tities related to the linear response to an applied volt- 
age bias V, i.e. the conductance and the torkance. In 
particular, it has been pointed out that, in deriving the 
spin current, the charge conservation through the system 
must be monitored in order to guarantee the conserva- 
tion laws. Indeed, the spin current is not conserved due 
to the presence of a spin-transfer torque acting on the lo- 
cal magnetization of the free- layer and thus a violation of 
the charge conservation law may artificially change the 
spin current gradient. The above change modifies the 
spin torque in a quantitative way. As for the observ- 
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FIG. 11: Parallel (dashed line) and perpendicular (full line) 
component of the spin torque Tx.n as a function of the in- 
terface potential z (zi = Z2 = zz = z) computed setting the 
model parameters as follows: e/A = 0.01, hz = 0.5, F = 0.5, 
T) = 1/200, e = 7r/2, kFd2 = 300, while fefrfi = 300 for the 
upper panel, kpdi — 380 for the middle panel and kpdi = 400 
for the lower panel. 



ables, in the second part of the paper, we derived the 
conductance and the magnetoresistance of a supercon- 
ducting spin- valve and analyzed all the relevant quantum 
size and coherent effects. Our analysis showed evidence 
of Andreev reflections in the subgap transport at varying 
the Zeeman interaction revealing the importance of the 
spin-flip processes. As for the magnetoresistance we an- 
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alyzed quantum size effects due to the superconducting 
layer thickness and showed that it displays a strong oscil- 
latory and non-monotonous behavior as a function of the 
interlayer width. A peculiar interplay between supercon- 
ducting and spin polarized transport properties becomes 
evident for thickness of the order of the superconducting 
coherence length. As a probe of the spin-polarized trans- 
port we analyzed the spin-torque in the linear response 
regime and characterized its behavior as a function of the 
interface transparencies and direction of the magnetiza- 
tion between the fixed and the free layer. It has been 
found that the torque and magnetoresistance are both 
strongly enhanced by a non-zero barrier height at the in- 
terfaces. Our analysis can provide an efficient way of de- 
tecting spin polarized transport in experiments on mag- 
netic/superconducting heterostructures helping some ba- 
sic understanding and stimulating further studies. 



Appendix A: Tensor product 

In this work the sign ^ is employed to define the Kro- 
necker product or tensor product of matrices. Given the 
matrices A and B the matrix C = A(E) B is obtained as 
follows: 



C = 



AiiB ... AinB 



ArnlB 



■^mn'^ 



(Al) 



where the size of ^ is ?7i x n. According to the above 
definition, provided that |e) = (1,0)* and |— ) = (0,1)*, 
we get, for instance, |e) ® |— ) — (0, 1, 0, 0)*. 



Appendix B: Scattering fleld in 
momentuni- representation 

Within the scattering approach one assumes that far 
from the scattering center the particle is free and its lin- 
ear momentum p = hk is a, good quantum number for 
labeling the scattering states. According to this, the 
scattering field can be expanded in the eigenstates of the 
linear momentum operator V = 'Y^a Ppp ^ rii3lsp{—ihdx)- 
The eigenstates of V in our tensor product notation are 
defined by 



V\^f}ak{x)) =hk\^p,kix)), 



(Bl) 



where I'ipakix)) = (V2^)"^|/3) ® |(T)e"w'="= and 1/^2^ 
is a normalization factor. These set of states satisfy the 
completeness relation: 

Y, J dk\-^f,,k{x)){^p.k{x')\ ^l4x4S{x ~ x'), (B2) 

Per 

where the identity operator is written as I4x4 = 
^p„Ppj3 ® |o')(o'|. The generic wave-function |^(x)) = 



TliB (T '^0o'(a;)|/3) ®\cr), can be written in the basis set of 
the eigenstates of V as: 

|*(x)) ^Y.1 ^^^Akm ® \a)e^^^^^, (B3) 

where the coefficients $/3cr(fc) = J dx'(l)j3cr{x')e~'^^f^^ , are 
related to the projection of |^(a;)) on the eigenvectors of 
V. 



Appendix C: Self-consistent determination of the 
chemical potential and the conductance tensor 

As described in the main text, in the non-symmetric 
case the chemical potential of the scattering region fig ^ 
(^1 -|- /i2)/2 — fi and thus a self-consistent computation 
of /is is required. Its calculation follows from the charge 
current conservation^^, J2i "^oi^^ f^s{V)) = 0. Since in 
principle, such a condition implies the solution of an in- 
tegral equation, a great simplification follows in the lin- 
ear response regime in the applied voltage bias V. In 
this case the charge current flowing through the «-th lead 
is obtained as li = J^jdijif^j ~ /^s), where gij is the 
conductance tensor and the charge conservation implies 
J2i h — 0. Solving the latter equation (Kirchhoff 's law) 
with respect to ^s we have: 



IJ-s 



L^ij 9ij 



(CI) 



Prom the equation above it immediately follows that in 
the case of a two-terminal2^ symmetric system (where 
5ii = 522) the chemical potential /x^ is bias independent, 
Us = M- More generically, one can analyze the potential 
drop to the left and right junction, i.e.: 



^ij - fis = qeVXj, 
where the coefficients A,- are function of 



below: 



9v 



Al =- 



512 + .922 

l^ij 9ij 



(C2) 
as shown 



(C3) 



■IV ■ 

921 +9ii 



Z-^ij 9ij 

Observing that A2 — Ai = 1 one correctly recovers that 
M2 — Ml = 9eV^(A2 — Al) = QeV, whilc for a symmetric 
system Xj = {—y/2. From the definitions above one 
immediately infers that the electrochemical potential of 
the scattering region /ig is displaced from n = (/ii -|-/i2)/2 
according to the relation: 



Us = IJ- + 



2 L Y.^^9^3 



(C4) 



Noticing that in the linear response regime It = GiV 
where Gi = X^i 9ij^j ^'^d using the expression above for 
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A, one obtains^^: 



G 



.922,911 ^ 921912 



l^ij 9ij 



(C5) 



For symmetric systems the relation above can be simpli- 
fied as 



Gs 



(511 -5i2)/2 



(C6) 



Let us note that Eg. dCGp and Eq. p^ in the main text re- 
produces the result given in Eq.(ll) of Ref . |24| using the 
Lambert's method. However, since we are considering a 
one-dimensional structure in place of the bidimensional 
one, the angular integration J dO cos(6')[- • •] is not present 
in our result. 



Appendix D: Boundary conditions of the scattering 
problem 

To determine the scattering matrix coefficients one has 
to use the mode-matching technique as formulated in 
the theory of quantum wave-guides. According to this 
method, the BdG equation is solved in each branch and 
the resulting eigenmodes are used to expand the scat- 
tering wave-function. Each wave-function is than deter- 
mined by imposing proper boundary conditions^-. E.g. 
in the presence of a single particle magnetic potential 



U{x) = j5{x)h ■ a, n = {nx,ny, Uz) being the unit vector 
describing the magnetization direction (|np = 1), the 
BdG wavefunction \E'(a;) — {uf{x),ui{x),v^{x),vi{x)Y 
must satisfy the following boundary conditions: 



*(a; = 0+) = ^{x = 0-) 



(Dl) 



2m"f 



d^-i^ix = 0+) - d:,^{x = 0-) = ^-4*(x = 0+), 



where the 4x4 matrix A is defined as follows: 
fi ■ a 



A = 



n-a* 



(D2) 



In the case of a non-magnetic potential U{x) = ^S{x)Isp 
the previous boundary conditions must be modified sub- 
stituting A with the 4x4 identity I4X4, i.e. ^— 7>l4x4- 
In the absence of potential, i.e. 7 = 0, the boundary 
conditions imply the continuity of the BdG wavefunction 
and its derivative. 
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